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We present general numerical methods to construct witness operators for entanglement detection 
j^jy and estimation of the fidelity. Our methods are applied to detecting entanglement in the vicinity of 

r—{ | a six-qubit Dicke state with three excitations and also to further entangled symmetric states. All 

our witnesses are designed to keep the measurement effort small. We present also general results on 
the efficient local decomposition of permutationally invariant operators, which makes it possible to 
00 ' measure projectors to symmetric states efficiently. 
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I. INTRODUCTION 

> ■ 

Ch ' n, 

Entanglement plays a central role in quantum mechanics and in quantum information processing applications 

Moreover, it appears also as the main goal in today's quantum physics experiments aiming to create various quantum 
states [2] . For example, entanglement has been realized with photonic systems using parametric down-conversion and 
conditional detection H, 3, H,H, Bi H, Hi with trapped cold ions [1(1 EH, G3 , in cold atomic ensembles in cold 
, atoms in optical lattices |14j and in diamond between the electron and nuclear spins [l5| . These experiments aimed at 
J> ■ creating entangled states. Entanglement makes it possible for some quantum algorithms (e.g., prime factoring, search 
in a database) to outperform their classical counterparts. Entangled particles are needed for quantum teleportation 
■ and other quantum communication protocols. Moreover, the creation of large entangled states might lead to new 
' insights about how a classical macro-world emerges from a quantum micro-world. 

In a multi-qubit experiment, typically the full density matrix is not known, and only few measurements can be 
made, yet one would still like to ensure that the prepared state is entangled. One possibility is applying entanglement 
witnesses [ill [I?} • These are observables that have a positive expectation value for separable states, while for some 
entangled states their expectation value is negative. Since these witness operators are multi-qubit operators, they 
typically cannot be measured directly and must be decomposed into the sum of locally measurable operators, which 
are just products of single-qubit operators pi [l8l Il9j . 

For many quantum states, like the Greenberger-Horne-Zeilinger (GHZ, [2(|) states and the cluster states [2l| such a 
decomposition of projector-based witness operators seems to be very difficult: The number of terms in a decomposition 
to a sum of products of Pauli matrices increases rapidly with the number of qubits. However, practically useful 
entanglement witnesses with two measurement settings can be constructed for such states [5|, [22J . It also turned out 
that there are decompositions of the projector for GHZ and W-states in which the increase with the number of qubits 
is linear [23| . 

However, optimal decomposition of an operator is a very difficult, unsolved problem. Moreover, in general, it is still 
a difficult task to construct efficient entanglement witnesses for a given quantum state. For that, typically we need 
to obtain the maximum of some operators for product states. In most of the cases we would like to detect genuine 
multipartite entanglement. For that, we need to obtain the maximum of these operators for biseparable states, which 
is again a very hard problem. 

In this paper our goal is to design witnesses that make it possible to detect genuine multipartite entanglement with 
few measurements, and also to estimate the fidelity of an experimentally prepared state with respect to the target 
state. Here three strategies are applied to find an experimentally realizable witness, (i) The first strategy is based on 
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measuring the projector-based witness 

W (P) = const. • 1 - |*)(*| (1) 

for the detection of genuine multipartite entanglement. \^} is the target state of the experiment. For reducing 
experimental effort, the aim is to find an efficient decomposition of the projector, (ii) The second strategy is to find 
a witness that needs much fewer measurements than the projector witness, but the price for that might be a lower 
robustness against noise. The search for such a witness can be simplified if we look for a witness W such that 

W - aW (p) > (2) 

for some a > 0. Such a witness is guaranteed to detect genuine multi-qubit entanglement. The advantage of this 
approach is that the expectation value of W can be used to find a lower bound on the fidelity, (iii) The third strategy 
is to find a witness independent from the projector witness. In this case one has to find an easily measurable operator 
whose expectation value takes its maximum for the target state. Then, one has to find the maximum of this operator 
for biseparable states. Any state that has an operator expectation value larger than that is genuine multipartite 
entangled. 

For the optimization of entanglement witnesses for small experimental effort and large robustness to noise, we use 
semidefinite programming [13, HE, [H, H3, HH Ull ■ Our methods can efficiently be used for multi-qubit systems with 
up to about 10 qubits. This is important, since there are many situations where semidefinite programming could help 
theoretically, but in practice the calculations cannot be carried out even for systems of modest size. 

We use our methods to design witnesses detecting entanglement in the vicinity of symmetric Dicke states. An 
./V-qubit symmetric Dicke state with m excitations is defined as [30j, l3l| 

^ ' k 

where ^uVki-) denotes summation over all distinct permutations of the spins. ') is the well known iV-qubit W 
state. The witnesses we will introduce in the following have already been used in the photonic experiment described 
in [32|], aiming to observe a \D^} state [H, l33| . We show that genuine multi-qubit entanglement can be detected 
and the fidelity with respect to the above highly entangled state can efficiently be estimated with two and three 
measurement settings, respectively. As a byproduct, we will also derive an upper bound for the number of settings 
needed to measure any permutationally invariant operator. We show that such operators can efficiently be measured 
even for large systems. 

The structure of our paper is as follows. In section 2, we present the basic methods for constructing witnesses. 
In section 3, we use these methods for constructing witnesses to detect entanglement in the vicinity of a six-qubit 
symmetric Dicke state with three excitations. In section 4, we present witnesses for states obtained from the above 
state by measuring some of the qubits. In Appendix A, we summarize the tasks that can be solved by semidefinite 
programming, when looking for suitable entanglement witnesses. In App endix B, we summarize some of the relevant 
numerical routines of the QUBIT4MATLAB 3.0 program package [3J|. In Appendix C, we present entanglement 
conditions for systems with 5 — 10 qubits that will be relevant in future experiments. 



II. BASIC DEFINITIONS AND GENERAL METHODS 



A multi-qubit quantum state is entangled if it cannot be written as a convex combination of product states. 
However, in a multi-qubit experiment we would like to detect genuine multi-qubit entanglement (3f|: The presence of 
such entanglement indicates that all the qubits are entangled with each other, not only some of them. We will now 
need the following definitions: 

Definition 1. A pure multi-qubit quantum state is called biseparable if it can be written as the tensor product 
of two, possibly entangled, multi-qubit states 

I*) = O |*2). (4) 

A mixed state is called biseparable, if it can be obtained by mixing pure biseparable states. If a state is not biseparable 
then it is called genuine multi-partite entangled. In this paper we will consider witness operators that detect 
genuine multipartite entanglement. 

Definition 2. While an entanglement witness is an observable, typically it cannot be measured directly. This 
is because in most experiments only local measurements are possible. At each qubit k we are able to measure a 
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single-qubit operator Mf., which we can do simultaneously at all the qubits. If we repeat such measurements, then 
we obtain the expectation values of 2 N — 1 multi-qubit operators. For example, for N = 3 these are Mi ® 1 ® 1, 1 ® 
M 2 ® 1, 1 ® 1 ® M 3 , Mi ® M 2 ® 1, Mi ® 1 ® M 3 , 1 ® M 2 ® M 3 , Mi ® M 2 ® M 3 . The set of single-qubit operators 
measured is called measurement setting [4| and it can be given as {Mi, M 2 , M 3 , M/v}. When we consider an 
entanglement condition, it is important to know, how many measurement settings are needed for its evaluation. 

Definition 3. Many experiments aim at preparing some, typically pure quantum state g. An entanglement 
witness is then designed to detect the entanglement of this state. However, in real experiments such a state is never 
produced perfectly, and the realized state is mixed with noise as given by the following formula 

£?noisy ^noisc 

) = (1- (5) 

where p no i S0 is the ratio of noise and g no i se is the noise. If we consider white noise then g no i SC — TL/2 N . The noise 
tolerance of a witness W is characterized by the largest p no isc for which we still have Tr(Wgnoisy) < 0. 

In this paper, we will consider three possibilities for detecting genuine multi-qubit entanglement, explained in the 
following subsections. Later, we will use these ideas to construct various entanglement witnesses. 



A. Projector witness 



A witness detecting genuine multi-qubit entanglement in the vicinity of a pure state | VP) can be constructed with 
the projector as 

W^ P) :=\%l-\*){% (6) 

where A is the maximum of the Schmidt coefficients for \ty), when all bipartitions are considered For the states 
considered in this paper, projector-based witnesses are given by [H, H3 , Wi\ 

^ V D(N,N/2) • 2iV-l JL \ U N I \ N I) \< ) 

<li) ^l-I^X^I- (8) 

These witnesses must be decomposed into the sum of locally measurable terms. For this decomposition, the following 
observations will turn out to be very important. 

Observation 1. A permutationally invariant operator A can always be decomposed as [45| 



A = J2c n aT, (9) 

n 

where a n are single qubit operators, and such a decomposition can straightforwardly be obtained. 
Proof. Any permutationally invariant multi-qubit operator A can be decomposed as 

A = H c ™ H P k( B n,l ® B na ® B nt3 ® ... ® B ntN )P k , (10) 
n k 

where B n ,m are single qubit operators, c„ are constants, and Pk are the full set of operators permuting the qubits. 
For odd TV, we can use the identity 

^P fc (B„,i ® B n .2 ® B„ j3 ® ... ® B 1hN )P k 
k 

= 2 -(N-i) ( SlBnA + S2 B n , 2 + S3 B rh3 + ...)® N . 

Si,S 2 , ... = ±1, 
S1S2S3 ■ ■ ■ SN = +1 

(11) 

Substituting (llip into (|10p , we obtain a decomposition of the form §9§ . Equation (jlip can be proved by carrying out 
the summation and expanding the brackets. Due to the sis 2 s 3 • ■ ■ sn = +1 condition, the coefficient of B n> i ® i? n>2 ® 
B n 3 ® ... ® B n jy is 1. The coefficient of terms like B n l ® B n 1 ® B n 3 ® ... ® B n ^, that is, terms containing one of 
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the variables more than once is zero. For even N, a similar proof can be carried out using [4£ 
J2 P k( B n,i ® B na ® B n<3 <g> ... <g> £„ 



= 2 -(iv-i) J2 s 1 (B nil + s 2 B ni2 + s 3 B r , 

Si, s 2 , ... = ±1, 
S1S2S3 • • • Sn = +1 

(12) 

Next, we give two examples for the application of (fTT|) and (fT2l for the decomposition of simple expressions 

^ftK^a^Pfc = ijc^ + a^-K-ov,)® 2 !, (13) 

Pk{e x ®a y ® a z )P k = \\ (a x + cr y + <t z )® 3 + (a x -cr y - a z ) m 

k ^ 

+ (~<J x ~a v +<J z )^ + (-a x +a y -a z f 3 }, (14) 



where are the Pauli spin matrices. While the first example does not reduce the number of settings needed, the 
second example reduces the number of settings from 6 to 4. 

Next, we present a method to get efficient decompositions for permutationally invariant operators. 
Observation 2. Any iV-qubit permutationally invariant operator A can be measured with at most 

C N = |iV 3 + iV 2 + fiV (15) 

local measurement settings, using (jlll) and ([12]) . 

Proof. We have to decompose first A into the sum of Pauli group elements as 

A= c ijm Y,Pk{aT®(7f®crf m ®l^ N - i -i- m: ))P k , (16) 

z,j,m: i-\-j-\-m<N k 

where Cij m are some constants. Then, such a decomposition can be transformed into another one of the form 
using (fTTj) and (fT2|) . All of the settings needed are of the form {a, a, a, a} where a = n x a x + n y a y + n z a z , rik are 
integer and 1 <J2 k \ nk \ — Simple counting leads to an upper bound Cn for the number of settings given in (fT5|) . 
Here we considered that (n x ,n y ,n z ) and (— n x , — n y , — n z ) describe the same setting. An even better bound can be 
obtained using that (n Xl n y ,n z ) and (cn x , cn y , cn z ) for some c ^ represent the same setting. An algorithm based 
on this leads to the bounds C' N = 9, 25, 49, 97, 145, 241, 337, 481, 625 for N = 2, 3, 10 qubits, respectively. 

For the projector |Z)jy^ 2 ' > )(Z)jy^ 2 ' > |, the decomposition to Pauli group elements contain only terms in which each 
Pauli matrix appears an even number of times. Hence, all of the settings needed are of the form {a, a, a, a} where 
a = 2n x a x + 2n y a y + 2n z a z , rik are integer and 1 < J2k \ n k\ ^ N/2. For this reason, £jv/2 an d £jv/2 are u PP er bounds 
for the number of settings needed to measure this operator. 

Let us discuss the consequences of Observations 1 and 2. They essentially state that the number of settings needed 
to measure a permutationally invariant operator scales only polynomially with the number of qubits. This is important 
since for operators that are not permutationally invariant, the scaling is known to be exponential (36| . Moreover, even 
if we can measure only correlation terms of the form o® , we can measure any permutationally invariant operator. 



B. Witnesses based on the projector witness 



We can construct witnesses that are easier to measure than the projector witness, but they are still based on the 
projector witness. We use the idea mentioned in the introduction. If is the projector witness and is fulfilled 

for some a > 0, then W is also a witness. This is because W has a negative expectation value only for states for 
which W( p ) also has a negative expectation value. The advantage of obtaining witnesses this way is that we can have 
a lower bound on the fidelity from the expectation value of the witness as 



Tr(e|*)(vl/|) > A 2 -iTr(W^). 



(17) 
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We will look for such witnesses numerically, such that the noise tolerance of the witness be the largest possible. 
This search can be simplified by the following observation. 

Observation 3. Since we would like to construct a witness detecting genuine multiqubit entanglement in the vicinity 
of a permutationally invariant state, it is enough to consider witness operators that are also permutationally invariant. 
Proof. Let us consider a witness operator that detects entanglement in the vicinity of a permutationally invariant 
state g and its expectation value takes its minimum for g. Then, based on ((5|) , the witness W detects entanglement if 

Tr(Wg) 

p ^> Tl{W g)-T,(Wg noisc )- (18) 

For a permutationally invariant state g, we have g = ^2 k PkgPk, where Np is the number of different permutation 
operators We assume that the same holds also for gnoise- Let us define the permutationally invariant operator 
W' = J2k PkWPk- The operator W' is non-negative on all biseparable states since 

inf Tr(Wg) = ^ V inf Tr{WP k gPk) < inf Tr(W'g), (19) 

geB k 

where B is the set of biseparable states. Hence, W' is a witness detecting genuine multipartite entanglement. Since 
we have Ti(Wg) = Tr(W"g), and Tr(yVp noiso ) = Tr(W"0 noisc ), the robustness to noise of W is identical to that of W. 
Hence, it is sufficient to look for witnesses that are permutationally invariant. 

We will first consider measuring the { o~ x , o~ x , a x , o~ x , a x , o~ x } and {o~ y , <r y , cr y , cr y , cr y , o~ y } settings, where oi are the 
Pauli spin matrices. This we call the two-setting case. Then we will consider measuring also the {<j z , a z , o~ z , o~ z , a Zl o~ z } 
setting, which we call the three-setting case. Due to Observation 3, we consider only permutationally invariant 
witnesses. Such witnesses can be written as 

JV 

W(a , {am}) := a • 1 + £ J2 a ^J2 Vk ^" ® *® {N ~ n) ], (20) 

l=x,y,z n—1 k 

where the summation is over all distinct permutations, and ao and ai n are some constants. We will consider a simpler 
but equivalent formulation 

N 

W(c ,{c /n }) :=c -a+ H c '» j r. ( 21 ) 

where Co, ci n are the coefficients of the linear combination defining the witness and J; are the components of the total 
angular momentum given as 

i N 

k=l 

Here erf denotes a Pauli spin matrix acting on qubit (fc). 

Finally, if we consider detecting entanglement in the vicinity of \D^^ 2 ^) states, then further simplifications can be 
made. For this state and also for the completely mixed state all odd moments of Ji have a zero expectation value. 
For any witness of the form (f2Tj) . the maximum for biseparable states does not change if we flip the sign of ci n for all 
odd n. Hence, following from an argument similar to the one in Observation 3 concerning permutational symmetry, 
it is enough to consider only even powers of J; in our witnesses. 

C. Witnesses independent from the projector witness. 

In general, we can also design witnesses without any relation to the projector witness. We can use an easily 
measurable operator M to make a witness of the form 

W := cl - M, (23) 

where c is some constant. To make sure that ()23[) is a witness for genuine multipartite entanglement, i.e, (W) is 
positive on all biseparable states, we have to set c to 

c = max (M)i$\, (24) 
|*>ee 1 ' 
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where B is the set of biseparable states. The optimization needed for (|24]) can be done analytically. For example, 
for the \D 4 ) state a witness has been presented that detects genuine four-qubit entanglement by measuring second 
moments of angular momentum operators (37j . However, analytical calculations become exceedingly difficult as the 
number of qubits increases. 

The optimization can also be done numerically, but one cannot be sure that simple numerical optimization finds the 
global maximum. (See Appendix B for a reference to such a MATLAB program.) Semidcfinitc programming is known 
to find the global optimum, but the optimization task (|24j) cannot be solved directly by semidefinite programming. 
Instead of looking for the maximum for biseparable states, using semidefinite programming, we can look for the 
maximum for states that have a positive partial transpose [HI, l3a |. (See Appendices A and B.) This way we can 
obtain 

c' := max max (M) e , (25) 

I g>0,g T i>0 

for which c' > c. The first maximization is over all bipartitions /. Thus, when putting c' into the place of c in (|23l) . 
we obtain a witness that detects only genuine multipartite entanglement. In many cases simple numerics show that 
c = c' . In this case our witnesses are optimal in the sense that some biseparable state gives a zero expectation value 
for these witnesses. 

Finally, let us discuss how to find the operator M in (|23[) for a two- or a three-setting witness, in particular, for 

detecting entanglement in the vicinity of Based on section lll Bl we have to look for an operator that contains 

only even powers of J/ . Hence, the general form of a two-setting witness with moments up to second order is 

Wd ( n,n/2) : = c dn - (J| + 4), (26) 

where cdn is a constant [39j]. The coefficients of and Jy could still be different, however, this would not lead to 
witnesses with a better robustness to noise. 

For other symmetric Dicke states, based on similar arguments, a general form of a witness containing moments of 
J i up to second order such that it takes its minimum for \D ™') is of the form 

W D(N, m) : = c i ~ (J" + Jy) + Q(J* - (^>|<->>) 2 ' ( 2? ) 

where c q and q are constants. For the witnesses described in this section, the optimization process is more time- 
consuming than for the witnesses related to the projector witness. Because of that we presented witnesses of the 
above type that are constructed only with the first and second moments of the angular momentum operators, and 
thus contain few free parameters. 

III. WITNESSES FOR A SIX-QUBIT DICKE STATE WITH THREE EXCITATIONS 

In this section, we will consider entanglement detection close to a six-qubit symmetric Dicke state with three 
excitations, denoted as \D^). There are several proposals for creating Dicke states in various physical systems 

0,111113,111. 

A. Witnesses based on the projector witness 

1. Two-setting witness 

Let us consider the two-setting case and define first the optimization problem we want to solve. We would like to 
look for the witness W with the largest noise tolerance that fulfills the following requirements: 

(i) W is a linear combination of certain basis operators Bk, that is, W = J2k c kBk, 

(ii) W - aW^l 3) > with some a > 0. 

For the two-setting case we set = {1, Jj, Jy, J^, Jy, J®, Jy}- The second condition makes sure that W is also a 

witness detecting genuine multipartite entanglement. 

Note that any optimization algorithm can be used for looking for W. Even if we do not find the global optimum, that 
is, the witness with the largest possible robustness to white noise, W is still a witness detecting genuine multipartite 
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0.2 0.4 D 0.6 0.8 1 

^noise 



FIG. 1: The fidelity F vs. noise (solid) and the fidelity estimate F' vs. noise (dotted), for the white noise (bottom two curves) 
and for the non- white noise (|30[) (top two curves). For the fidelity estimate, the three-setting witness ((29} was used. 

entanglement. However, semidefinte programming can be used to find the global optimum. (See Appendix A.) The 
two-setting witness obtained this way is 

Wggg := 7.75 ■ 1 - f§ (J 2 X + 4) + § (J x 4 + Jy) - ^(J e x + J e y ), (28) 

which tolerates white noise if p no i SC < 0.1391. Straightforward calculation shows that W^g gs — 2.5W ( ' P ' ) > 0. Based 
on HO), (0.6 - (W^ 6 2) 3) )/2.5) bounds the fidelity from below. 

2. Three-setting witness 

Similarly we can look for the optimal witness for the three-setting case. The result is 

^D(6,3) : ~ ' ^ ~~ TE^x Jy) ^ei^x + Jy) ~ TEo(Jx Jy) 

, 1007 t2 _ 31 t4 i _23_ 76 /oq\ 
' 360 J * 36 J z 360 J z- 

White noise is tolerated if p n oise < 0.2735. It is easy to check that W is a witness as W^g L — 2.5W' P ' > 0. 

Based on (|T7|). the expectation value of this witness can be used to bound the fidelity as F > 0.6 — (VV^g 3 \)/2.5 =: 
F' . Here we will demonstrate how well the fidelity estimation works for our witness for noisy states. We consider first 
white noise, then non-white noise of the form 

ffi? i=PD*\D?){D?\ + i=£» (lufWl + l4 4) )(4 4) l) , 00) 

with PD63 = 4/7, which is one of the relevant types of noise for the experiment of (32|- Note that the noise contains 

the original state \Dg'){D£'\. The results are shown in figure [1] For the non- white noise (|50|) . the fidelity estimate 
based on the witness yields a very good estimate. 

Note that it is also possible to design a witness for the largest possible tolerance to the noise in ([30|) . Due to the 
special form of the noise, the fidelity estimate turns out to be equal to the fidelity. This is remarkable: The fidelity 
can be obtained exactly with only three local measurements. 

3. Measuring the projector-based witness 

For measuring the projector-based witness ([7]) for N — 6, one has to decompose the projector in an efficient way. 
The straightforward decomposition into the weighted sum of products of Pauli spin matrices leads to a scheme that 
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(a) (b) 

FIG. 2: (a) The measurement settings needed to measure the projector to the six-qubit symmetric Dicke state with three 
excitations based on the decomposition (|3ip . A point at (x,y,z) indicates measuring xa x + ya y + za z on all qubits. (A) 
Settings for Mermuiz.z, (v) settings for Mermin H , z , (□) a x ± a y ± a z , and (O) res t of the settings, (b) Settings for the 
four-qubit Dicke state with two excitations based on (I37|l . (□) ±a\ T , ±<J y , ±ct z , and (Q) o- x ± a y ,a x ± a z , and a y ± a z . 



needs 183 settings, since for all local operators all the permutations have to be measured. The number of settings 
needed can dramatically be decreased if one is looking for a decomposition of the form Observation 1 makes it 
possible to decompose the projector in this way such that only 25 settings are needed. We could further decrease the 
number of settings needed and found the following decomposition 



(3)wr>(3)| 



M\D™){D t 



-0.6[1] + 0.3[x ± 1] - 0.6[a;] + 0.3[y ± 1] - 0.6[y] 
+0.2[z ± 1] - 0.2[z] + 0.2Mermino, z 
+0.05[x ± y ± 1] - 0.05[a; ± z ± 1] - 0.05[y ±z±l] 
-0.05[x ± y ± z] + 0.2[x ± z] + 0.2[y ± z] + 0.1[x ± y] 
+0.6Mermin K , + 0.6Mermin 



(31) 



Here we use the notation [x + y] 
over the two signs, i.e., [x ± y] = 



[x + y + 1] = (<j x + o y + l)® 6 , etc. The ± sign denotes a summation 



[x + y] + [x — y] . The Mermin operators are defined as 



Mermin a , 6 := ^ {-l) k ' 2 ^ T k (e. 



-1 a b), 



(32) 



k even 



where <tq = 1. That is, it is the sum of terms with even number of <t q 's and cr^'s, with the sign of the terms depending 
on the number of cr a 's. The expectation value of the operators Mermin a b can be measured based on the decomposition 

M 



N 



Mermin a ,h = 



N 



k=l 



cos 



(33) 



Hence, Mermin XjZ and Merminj,^ can be measured with six settings. Mermino !Z , on the other hand, needs only 
the measurement of the {a z , a Zl cr z , er z , a z , a z } setting. Knowing that [A], [A + 1] and [A — 1] can be measured 
with a single setting {A, A, A,..., A}, we find that 21 measurement settings are needed to measure \Dq )(Dq \ : 
x, y, z, x ±y,x ± z,y± z, \/3x ± z, \[Zz ± x, V3y ± z, ^/3z ± y, and x±y±z. The settings are also shown in figure [21(a) 



B. Witness independent from the projector witness 

So far we constructed witnesses that detected fewer states than the projector-based witness, in return, they were 
easier to measure. When proving that they were witnesses, we used the simple relation @. Following the example 
of [37[, we now look for a two-setting witness of the form ([26]) for N = 6 that is independent from the projector 
witness. For determining Coq, we need to compute the maximum of J| + Jy for biseparable states for all the possible 



9 



bipartitions. As we have discussed in section Hi CI instead of looking for the maximum for states that are separable 
with respect to a certain bipartition, we can also look for the maximum for PPT states. (See Appendix A.) We obtain 

c m := 11.0179. (34) 

(12) ... ... n . 

Ww 6 3 ) detects genuine multipartite entanglement if for white noise p n oisc < 0.1091. Simple numerical optimization 
leads to the same value for the maximum for biseparable states [5l| . Hence we find that our witness is optimal. 
Finally, the list of witnesses presented in this section are shown in the top part of table 1. 

TABLE I: The list of entanglement witnesses presented in this paper, together with the number of measurement settings needed 
to measure them and their robustness to white noise. Top four lines: six-qubit witnesses. Bottom five lines: four- and five-qubit 
witnesses. 



Witness Number of settings Noise tolerance 



<i,3) 


21 


0.4063 


<a 3) 3, 


3 


0.2735 


KV D(6,3) 


2 


0.1391 


KK D(6,3) 


2 


0.1091 




2 


0.1046 




7 


0.2667 


Wd ( 4,i)(?=1-47) 


3 


0.1476 


><i, 2 ) 


9 


0.3556 


w (P3) 


2 


0.2759 



IV. WITNESSES FOR STATES DERIVED FROM \D£>) VIA PROJECTIONS 

By projective measurements of one or two of the qubits we can obtain several states that are inequivalent under 
stochastic local operations and classical communication (SLOCC). Surprisingly, these states still possess genuine 
multipartite entanglement [H,[44]]. Next, we discuss how to detect the entanglement of these states. 

A. Witnesses for the superposition of five-qubit Dicke states: 

After measuring one of the qubits in some basis and post-selecting for one of the two outcomes, one can obtain 
states of the form 

g V5 := Cl \Di 2) )+c 2 \Di 3) ), (35) 

where |ci| 2 -I- | C2 1 2 = 1- For such states, the expectation value of + Jy is maximal, thus a witness of the form 
for N = 5 is used to detect their entanglement. Both semidefmite programing and simple numerical optimization 
leads to cds '■= 7.8723. Naturally, (W^g ) is minimal not only for states of the form (f35|) . but for any mixture of such 
states. 



B. Witness for the four-qubit W-state 

Now we will construct witnesses for a four-qubit W-state, which is obtained from \D§ ) if two qubits are measured 
in the a z basis, and the measurement result is +1 in both cases. We consider a witness of the form (|2"T)) for N = 4 and 
m = 1. We try several values for q and determine c q for the witness Wd(4,i)(<z) as a function of q using semidefmite 
programming. For each witness we also compute the noise tolerance. The results of these computations can be seen 
in figure [31 It turns out, that the best witness is obtained for q — 1.47 and c q — 4.1234. It tolerates white noise if 
p noiso < 0.1476. 
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C. Three-setting witness for the four-qubit Dicke state 

(2) (3) 

A \D\ ) state can also be obtained from \D G ), namely if the measurement outcomes are +1 and —1 for two 
consecutive u z measurements. For that case, we look for a three-setting witness, based on the projector witness. For 
white noise, the result is 

HgJ 2) := 2 • 1 + |(J X 2 + J 2 y - 4 - J y 4 ) + fl Jl - ^J A Z . (36) 
The witness tolerates white noise if p n oisc < 0.2759. It is easy to check that W is a witness: One has to notice that 

( P3) ( P) (P) XF^Y, 

^D(4 2) ~ 3^"d(4 2 ) > 0, were ^ is defined in ([7]). Thus, the fidelity can be estimated from the measurement of 

the witness as F > 2/9 - (W^ ( f 2) )/3. 



D. Measuring the projector witness for the four-qubit Dicke state 



We can also measure the projector witness W^j 2 ^ = |l — \D^) (D^ \ . The method in Observation 1 gives the 
following decomposition for the projector 

16|Z?f)<^ 2) | = l([x] + [x±l] + [y] + [y±l}) 

+ l(8[z}-[z±t]-[x±z]-{y±z]) + ±{x±y]. (37) 

The 9 measurement settings are x, y, z, x±y, x±z, and y ±z, shown also in figure[2Ib). The list of witnesses presented 
in this section are given in the bottom part of table 1. 



V. CONCLUSIONS 



In summary, we presented general methods for constructing entanglement witnesses for detecting genuine mul- 
tipartite entanglement in experiments. In particular, we considered projector-based witnesses and found efficient 
decompositions for them. Then, we constructed two- and three-setting witnesses for symmetric Dicke states that 
were based on the projector witness, as well as independent from the projector witness. We applied our methods to 
design witnesses for the recent experiment observing a six-qubit symmetric Dicke state with three excitations [32l ]. 
Our methods can be generalized for future experiments. As a first step, in Appendix C we list some entanglement 
witnesses for systems with 5 — 10 qubits. Moreover, recent results on the symmetric tensor rank problem suggest that 
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decompositions more efficient than the one in Observation 1 are possible, however, they involve complex algorithms 
[45l ] . Thus, it would be interesting to look for better upper bounds for the number of settings used for symmetric 
operators. 
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APPENDIX A: SEMIDEFINITE PROGRAMMING USED FOR OBTAINING WITNESSES 

Here we summarize two optimization problems that are useful for designing entanglement witnesses and can be 
solved with semidefinite programming. Both tasks are related to designing witnesses that are easy to measure. 

1. Semi-definite programming can be used to find the witness W with the largest noise tolerance as explained in 
the beginning of section IIII A II The corresponding task can be formulated as 

minimize J2 k c k^{Bk Pnoisc ) ) 

subject to J2 k ckTr(B k g) = -1, (Al) 
E fc c k B k - aWC) > 0, 
a > 0. 

Here g is the state around which we detect entanglement. g no i SC is the noise, not necessarily white. The optimization 
is over a and the c k 's. 

2. Semi-definite programming can be used to look for the maximum for PPT bipartite states. This gives an upper 
bound on the maximum for biseparable states. In many cases, the two coincide. The corresponding task can be 
formulated as a standard semidefinite program as 

minimize — Tr(Mg), 

subject to g > 0, (A2) 
Tr(g) = 1, 
Q Ta > 0. 

Here Ta means partial transpose according to some groups of the qubits. 

APPENDIX B: LIST OF MATLAB SUBROUTINES 

We summarize some of the MATLAB routines of the QUBIT4MATLAB 3.0 package that can be used for the 
calculations necessary for designing entanglement witnesses. A full list of the commands is given in [34| . 

The command decompose can be used to obtain a decomposition of a Hermitian operator into the sum of products 
of Pauli spin matrices. Moreover, maxsep and maxsymsep can be used for getting the maximum for separable multi- 
qudit states and symmetric product states for a Hermitian operator, respectively. The command maxbisep gives the 
maximum for states that are biseparable with respect to some partitioning of the qubits. The command maxb gives 
the maximum for all possible bipartitions. All these commands look for the maximum with a simple optimisation 
algorithm that is not guaranteed to find the global maximum, nevertheless, it typically does find it. overlapb gives 
the maximum overlap of a state \^/) and biseparable states. It can be used to construct entanglement witnesses of 
the type ©. 

For semidefinite programming, we used SeDuMi [46| and YALMIP [47|. Two subroutines based on them are now 
in QUBIT4MATLAB 4.0 (48[. The command optwitness looks for the best witness that can be composed linearly 
from a set of operators, while maxppt determines the maximum of an operator expectation value for states with a 
positive partial transpose for some bipartitioning of the qubits. 
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APPENDIX C: WITNESSES FOR SYSTEMS WITH 5-10 QUBITS 
A three-setting witness based on the projector witness for the state \D^} is given by 

4 

Wg 4) := 1.3652-1+ ]T $>™ J '"> (CI) 

l=x,y,z n—1 



(C2) 



with 

/ 0.0038612 -0.0052555 0.0015016 -0.00010726 ^ 

{cin} = 0.0038612 -0.0052555 0.0015016 -0.000107266 . 

\ 3.124 -1.07699 0.11916 -0.0038992 j 

The noise-tolerance for white noise is p no i S e < 0.2578. For larger A, we can use the ansatz 

N/2 

W%$,n/2) : = c i ' 1 + C *y{^* + l ) 9N + ^ ~ l ) 9N + K + l) t3W + iPv - l) 0Ar } + E c - J *"- ( C3 ) 

n=l 

For A = 10, the optimal coefficients are c\ = 1.3115, c xy = —0.0023069, and c z = 
{3.4681,-1.2624,0.16494,-0.0084574,0.000146551}. White noise is tolerated if p noise < 0.2404. The large 

noise-tolerance for the A = 10 case suggests that a robust three-setting witness for \D^^) might be constructed 
even for large A. 

A three-setting witness independent of the projector witness for the A-qubit VF-state is given by (|27|) for m = 1. 
For A = 5, we have C5 = 5.6242, q$ = 2.22, and the witness tolerates white noise if p n oisc < 0.0744. For A = 6, we 
have C6 = 7.1095, q$ = 3.13, and noise is tolerated if p no iso < 0.0401. 
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